Quark contribution to r-mode instabilities for several phases of deconfined quark 

matter 
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R-mode instabilities lead to specific signatures in the evolution of rotating pulsars, and may 
provide a unique signature which can identify the phase of matter that exists in the interiors of such 
objects. 

The contributions of quarks' bulk and shear viscosities to the dissipation of r-mode instabilities are 
studied. It is shown that (contrary to earlier works) the quark contribution is not fully suppressed 
in the presence of fully gapped color-superconducting phase due to the fact that the bulk viscosity 
is a resonance effect between the density oscillation on one side, and the interaction rates that try to 
restore /3-equilibrium during these oscillations on the other. This will have an effect on the structure 
of the r-mode instability window and will affect the angular momentum evolution of pulsars. 

PACS numbers: 12.38.-t, 12.38.Aw, 12.38.Mh, 26.60.+C 



I. INTRODUCTION 

Compact (neutron) stars provide a natural laboratory 
of matter under extreme conditions. In the central re- 
gions of such stars the baryon density of matter could 
reach values up to 10 times the nuclear saturation den- 
sity (i.e., lOpo where po — 0.15 fm~ 3 ). At such high den- 
sities matter is likely to be in a deconfined state in which 
quarks, rather than hadrons, are the natural dynamical 
degrees of freedom 

0, H H II H ■ At sufficiently low tem- 
peratures, the ground state of deconfined quark matter 
is a color superconductor. (For reviews on color super- 
conductivity, see Refs. [1, 01,1, [Mill El El)- Many 
phases of color superconductivity are known that could 
possibly be realized in dense matter. It remains unclear, 
however, which of these describe the ground state of mat- 
ter under the specific conditions in stars. This is because 
of theoretical uncertainties in treating the strongly cou- 
pled, non-perturbative dynamics in QCD at the baryon 
densities of relevance. 

There are many instabilities associated with a star's 
rotation that can be operative in a neutron star. The 
most familiar of those is that if the star rotates with 
frequencies above a critical value, known as the Kepler 
frequency, mass shedding sets in from the star's equator. 
This frequency is i>k ~ 1200Hz in a typical neutron star. 

Of concern to us are the r-mode instabilities. (For re- 
views on the r-mode instabilities see REfs. [Hill). It- 
modes are fluid waves in rotating stars with the Coriolis 
force acting as the restoring force. These waves couple to 
gravitational radiation emission causing the star to lose 
angular momentum. Only when dissipative phenomena 
dampen these r-modes can the star rotate without losing 
any angular momentum. This sets a limit on the maxi- 
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mum rotation frequency of the star at any given temper- 
ature. This also affects the rotation frequency evolution 
of the star, see, for instance Ref. (l6j . 

There can be many dissipative phenomena that will 
damp the r-modes, these phenomena depend strongly on 
the phase of matter expected and on the structure of the 
star. In this paper we study the quark contribution to 
the r-mode instability window of quark stars for different 
phases of deconfined quark matter. 

At first, we consider the r-mode instability window for 
unpaired quark matter in comparison with the r-mode in- 
stability window of proton- neutron-electron matter (npe 
matter), the results are very similar to the earlier work 
in Ref. 0. 

It was suggested [I3| that the gapped quark contribu- 
tion in the color-flavor locked (CFL) phase is strongly 
suppressed even if the energy gap is as small as Acfl ~ 
lMcV. This will be reconsidered more carefully in this 
article. Moreover, the r-mode instabilities for the two- 
flavor color superconductor (2SC) phase will be studied, 
the results are different from those of Ref. [l7[ due to the 
fact that the bulk viscosity may increase with the reduc- 
tion of weak interaction rates as was seen in Ref. [HI]. 
The quark contribution to the r-mode instability window 
in the color-spin-locked (CSL) phase will also be studied. 

This paper is organized as follows: A quick review of 
the r-mode instabilities is presented in Sec. [ill Quarks 
contribute to the r-mode instability window mainly 
through two dissipative phenomena, the bulk viscosity 
and the shear viscosity, these are reviewed in Sec. IIIII 
In order to relate the dissipative factors "the bulk and 
shear viscosities" to the r-mode instabilities, one has to 
calculate the dissipation timescales related to these fac- 
tors, these are discussed in Sec. IIVI The results for the 
r-mode instability window are shown in Sec. El The dis- 
cussion of the results is given in Sec. ED 
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II. R-MODES 



Non-radial pulsations of the star couple to gravita- 
tional radiation (GR) emission and the GR emitted car- 
ries away energy and angular momentum from the star. 
In a non-rotating star, the effect of emitted GR is to 
dampen these oscillations. In a rotating star, however, 
the situation is rather different. The emission of GR 
causes the modes to grow, the reason being that modes 
which propagate in the direction opposite to the star's 
rotation have negative angular momentum as seen in the 
co-rotating frame. Therefore, GR lowers the already neg- 
ative angular momentum, i.e., the mode accelerates in 
the opposite direction of the star's rotation and the en- 
ergy of the mode grows, i.e., the amplitude of the mode 
grows. 

R-modes are fluid waves on rotating stars with the 
Coriolis force acting as the restoring force. The r-modes 
are primarily velocity perturbations, which are solutions 
to the perturbed fluid equations. For a slowly rotating 
star they have the form 0, El El III 



5v = aRtt 



B itot 



(i) 



where R and £1 are the radius and angular velocity of the 
unperturbed star, a is an arbitrary dimensionless ampli- 
tude of the mode, and is the magnetic-type spherical 
harmonic defined by 



y b 



rV x (rVYi„ 



(2) 



Because the Coriolis force dominates, the frequencies 
of the r-modes are independent of the equation of state 
and are proportional to the angular velocity of the star 

BE El, Hi 



(m- 1)(to + 2)„ 
ui = S2. 

TO + 1 



(3) 



The expressions for Sv and u) are the lowest-order terms 
in Q. 

Since r-modes are generally velocity perturbations, the 
energy of the mode has both kinetic and gravitational 
potential energy terms. We follow the Cowling approx- 
imation, in which the perturbation of the gravitational 
potential can be neglected. This has been shown to 
be a very good approximation, see, for instance, Refs. 
0, EE El H3] • The energy of the mode (measured in 
the co-rotating frame) is 



E 



p(Sv* ■ 5v)d 3 x. 



(4) 



This energy is conserved in the absence of dissipation. 
Including dissipation it satisfies [13, EH, El 



dE 
~dt 



oj(u) + mCl) ^ Ni 
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\SD lm \ 2 + 



U\5J lm \ 2 
c 2 (Z + l) 



[2n(5a* ab 5a ab ) + ((5<r*6a)] d 3 x. 



(5) 



The second term in Eq. ([5]) represents the dissipation 
due to viscosities of the fluid. It can be seen from the 
negative sign in front of the term that the effect is always 
to reduce the energy of the mode converting it to heat. 
The thermodynamic functions 77 and £ are the shear and 
bulk viscosities of the fluid. The viscous forces are driven 
by the shear 5a a b and the expansion scalar 5a of the 
perturbation defined by 



5a = W a 5v a . 



(7) 



The first term in Eq. represents the effect of GR; 
5Di rn and 5 Ji m are the mass and current multipole mo- 
ments of the perturbation, 



5D lm = I 5pr l Y l * m d 3 x, 



6 Jim = / r\p5v + 5pv) ■ Y^d 6 x 



(8) 
(9) 



The constants Ni 



- 4?rG w (i+lXj+2) nnsitive 

n — c 2i + i (i(z_i))[(2Z + l)!!p dlL poorcive 

and the sum is also positive definite (Gjv = 6.6742 x 

I0~ n m 3 Kg -1 s -2 is the gravitational constant). Thus, 

the effect of GR is determined by the sign of lo(uo + mH,), 

which is the product of the frequencies of r-modes in the 

inertial and the rotating frame, 



2(to- 1)(to + 2) 
(to + I) 2 



n 2 < 0, 



(10) 



which implies that the total sign of the first term in Eq (JSJ) 
is always positive. GR always increases the energy of the 
mode. 

It has been shown [3, El, [H, Hl[ that only modes with 
I = to can exist and that the most unstable mode is the 
one with I = to = 2, which will be the one considered in 
our calculations. 



III. DISSIPATIVE PHENOMENA 

A. Bulk viscosity 

In general, the bulk viscosity is a measure of the ki- 
netic energy dissipation during expansion and compres- 
sion of a fluid. In compact stars, the density oscil- 
lations of interest have characteristic frequencies that 
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ji s 







FIG. 1: Diagrammatic representation of the weak processes 
that contribute to the bulk viscosity of quark matter in stellar 
cores. 



J 



are of the same order of magnitude as the stellar ro- 
tation frequencies. These are bound from below and 



from above, 1 s 



i < 



< io 3 



(For the fastest- 



spinning pulsar currently known, PSR J1748-2446ad, one 
has oj w 4.5 x 10 3 s _1 corresponding to v = 716 Hz 
[22|). Under conditions in stars, the bulk viscosity of 
quark matter is determined by the combined effect of 
the flavor-changing weak processes. When an (instanta- 
neous) departure from chemical equilibrium is induced 
by expansion/compression of matter (Sfii = fi s — /id) and 
5fi2 = A*s — M« — Me), the weak processes try to restore the 
equilibrium state (Sfii = 0) and, while doing this, reduce 
the oscillation energy. 

These weak interaction are: the non-leptonic weak 
processes u + d <-> u + s shown diagrammatically in 
Figs- HI a) andQJb), and the Urea (semi-leptonic) pro- 
cesses d(s) — » u + e + D e , u + e — » d(s) + v e , see Figs.QJc)- 

mf). 

For studies of the viscosity in various phases of dense 

JUL 



nuclear matter, see Refs. [1 
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I3fl, 131|, 133, |33J, |34J, |3a 136|, HMMHEP' 

The bulk viscosity due to weak interactions in quark 
matter is given by, see Ref. [4l[ for a detailed derivation, 

C = Ci + C 2 + C 3 , (ii) 

where 



Ci 

C2 

Cs 



n ctict3 

u 9i + 92 

n a.ia.2 

ug\+ 9l 



aia 2 a 3 Cl + (ct\ + a 2 + a 3 ) {A 1 C 2 - A 2 Ci)' 
u\oi 2 azC\ + (ai + a 2 + 013) [(A 2 - B 2 ) C\ - A 2 C 2 ] 2 
ol\ol 2 oiz (Ci — C2) 2 + (ai + a 2 + as) {B\C 2 - B 2 C\ 

I 



(12a) 
(12b) 
(12c) 



with a.i = nto/Xi, n is the baryon density, to is the oscil- 
lation frequency, and Ai is the difference of the forward 
and reveres weak interaction rates to first order in 5/ii. 



three pairs of weak processes in Fig. [T] 



The thermodynamic functions Ai, Bi, and Ci are 
function of the chemical potentials and their derivatives 
with respect to the baryon number density, they are 
discussed in detail in Ref. 41], they do not vary sig- 
nificantly whether in the normal phase or in a color- 
superconducting phase as was seen in Ref. [IH . 



T(a) 


-r (6) = 


-XiSfii, 


r (c) 


-r (d) = 


-\ 2 s^ 2 , 


r (e) 


~ r (/) = 


-A 3 {5^ 2 - 6 hi) 



(13a) 
(13b) 
(13c) 



In order to further proceed with the calculation of the 
bulk viscosity in the normal phase of three-flavor quark 
matter, we also need to know the rate difference of the 



The rates of both the non-leptonic and semi-leptonic 
processes have been calculated in the literature, see 
Refs. jH S3, S3, SI and Refs. [HSl, respectively. In 
the limit of three massless quarks, for example, the rates 
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arc 



Ai 
A 2 
A 3 



P=Gtj, cos 2 9 C sin 2 e c ^T 2 , 

1 • r<1 ■ 2 a 2rri 

— — G F sm O c ^ s m s T , 

407T 



17 
15^ 



G F cos 2 9 c a s p d PuPeT 4 



(14a) 
(14b) 
(14c) 



Before we proceed, it is instructive to discuss the 
effect of color superconductivity on Aj's. Recall that 
Ti^SiH, A) - TiiS/H, A)_ = Xi(A)5m, see Eqs. (SI 
and that r^J/i^, 0) = r<(— tfytj, 0) where Ti(Sfj,i, A) 
is the forward interaction rate as a function of both 
(5/ii and A. Ti(5p,i, A) is the inverse interaction rate 
(backwards channel). For instance, if we choose IV a ) in 
Fig. [T] and Eqs. (fT3|) as the forward channel Fi(8fJ,i, A), 
Ti would be Fry in Fig. Q] and so forth. From the 
two equations mentioned above it is easy to show that 

r<(fyi,o) = -fi(6m,o) = |Ai(o)5 M i- 

Color superconductivity introduces a reduction factor 
Hi(A) to the interaction rate Ti and another suppression 
factor Hi(A) to the inverse interaction, 



TiiSnuA) = ri(8iii,0)Hi(A), (15) 
f t (SiM,A) = Ti(8iM,0)Si(^). (16) 

The relation between the gapped Aj(A) and the un- 
gapped Aj(0) becomes, 

A 4 (A)<5 M = Ti(SiM, O)ffi(A) -Ti(6iH, 0)H % (A) 
= Ti{Sin, 0)H l (A)-T l (~d^ O)ffi(A) 



g ,(A)+a(A ) <)t 



The effect of color superconductivity is to suppress the 
quark excitations exponentially, if all quark modes are 
gapped, we can neglect the quark shear viscosity from 
consideration since the shear viscosity will be reduced by 
approximately (including screening) e~sr. In this case, 
the quark contribution to the r-mode instabilities is solely 
governed by the bulk viscosity. Otherwise, only ugapped 
quarks' contribution will be considered. 

For the remainder of this article the values a s = 0.1, 
n = 5 pa will be used. 



IV. DISSIPATIVE TIMESCALES 

In order to compare the relative strengths of dissipative 
phenomena and GR emission it is convenient to calculate 
the dissipative timescales associated with these effects 
defined as 



2E 



(dE/dt)i 



(20) 



where the index i can stand for any dissipative phe- 
nomenon such as bulk viscosity (bv), shear viscosity (sv), 
or GR emission (GR). 

The energy of the mode can be found u sing Eqs. ([1]) 
and ([!]). Assuming spherical symmetry [bH. Il5l. Tl9j. 



E = -a 2 n 2 R- 2l + 2 
2 



/ pr 21 ^ 



2 dr. 



(21) 



(17) 



Gravitational Radiation Timescale 



Then the reduction to (0) would be a function of the 
suppressions to both the forward and the inverse inter- 
actions. 



Ai(A) = Ai(0) 



Hi(A) + Hi(A) 



(18) 



In the following, we will approximate the reduction 
H(A) to be of an exponential form for each incoming 
gapped quark branch in the interaction. 



B. Shear Viscosity 

The shear viscosity of ungapped quark matter is dom- 
inated by quark-quark scattering, for ungapped quark 
matter it is given by [H |47| 



i] = 5 x 10 



15 



0.1 



14/9 , rp \ -5/3 

(Wk) 



gem s . 



(19) 



The lowest-order contribution to the GR terms in 
Eq. ([5]) comes entirely from the current multipole mo- 
ment SJu. This can be evaluated using Eqs. ([TJ) and 
to lowest order in ft 0, El El 



SJu 



2afl / I 



cR 1 - 1 V I + 1 



pr 2l+2 dr. (22) 



This leads to a dissipative timescale [3, ES HH HI! 
1 32^Gf> 2 <+ 2 (l-l) 2 l fl + 2^ 2l+2 ' R 



T~GR 



c 2l+3 [(2/+i)!!]2 + i 



pr 



2l+2 dr. 
(23) 



Note that we are only interested in the case where I = 
m = 2. Assuming a uniform density star the dissipative 
"growth" time for GR emission becomes 



1 / CI 



1 



t gr 3.26 V irGp 



(24) 
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B. Shear Viscosity Timescale 



The time derivative of the mode energy due to the 
shear viscosity coefficient can be calculated from Eqs. ((T|), 
©, and 0, which leads to [II (US [HI 



(I-1)(2I + 1) . 
Using / = 2 and constant density, wc get, 



(25) 



1 dE 
E dt 



1 

TGR 



E 



Tdi 



(31) 



Where — 

Td 



— X— is the sum of the inverse of the dis- 
sipative timescales (but not the gravitational radiation 
timescale). Above the critical frequency the growth of 
the r-modes due to GR emission is the most dominant 
effect and the star emits GR and loses angular momen- 
tum until the effect of viscosities is the dominant one. 



1 _ 28tt t]R 
where M is the mass of the star. 



(26) 



C. Bulk Viscosity Timescale 

The time derivative of the co-rotating frame energy E 
due to the effect of bulk viscosity is 



(27) 




— =-U\V-8Wd?x. 



The reduction of dE/dt to a one-dimensional integral 
is not straight forward; the expansion scalar of the mode 
V • Sv is a complicated function of radius and angle. How- 
ever, since the bulk viscosity generally has no angular de- 
pendence, we may convert Eq. (|27p to a one-dimensional 
integral by defining an angle averaged expansion scalar 

squared ( V • Sv\ 2 




= -4tt I C (\^ ■ H ) r d r - (28) 



The angle-averaged scalar expansion is in general a 
complicated function, it has only been determined nu- 
merically in Ref. [HI, however, the simple analytical ex- 
pression 



V • Sv\ : 



a 2 fl 2 / r \ 6 
690 \RJ 



l + 0.86(- 



Q 2 \ 2 
ttGnpJ 



is an excellent fit to those numerical solutions [31| 
Using Eqs. (EH), (gg), and (EH) we get 



J_ _ {dE/dt) bv 

'Tbv E 



4.828 x 10" 



n 2 



ttGnP J pR 2 



(30) 

The critical angular velocity above which GR has the 
shortest timescale is found by solving [3, [H, [lj| 



V. THE R-MODE INSTABILITY WINDOW FOR 
QUARK STARS 

The dominant dissipative phenomena in a quark star 
are the quark bulk and shear viscosities. There are other 
dissipative phenomena that can play a role in determin- 
ing the r-mode instability window, for example, the shear 
viscosity due to electron-electron scattering, or, in the 
case of a superfluid phase, the friction between the nor- 
mal and the superfluid component that adds other vis- 
cosity coefficients [4{|. We are concerned with the con- 
tribution of quarks to the r-mode instability window. 

Our results are also of relevance to hybrid stars. How- 
ever, to get a complete picture, the effect of hadrons need 
to be considered as well, and friction between the quark 
core and the hadronic mantle will also play a decisive role 
in finding the r-mode instability window. 



R-mode Instability Window for the Normal 
Phase 



Using Eqs. and (0 alongside with Eqs. (ED, 1(21)]). 
and (|30[) one can calculate the R-mode instability window 
from Eq. (|3Tj) . The results for the critical frequencies are 
shown in Fig. [21 the region above the curve is the unstable 
region towards the emission of GR from r- modes. For 
comparison, the r-mode instability window for neutron- 
proton-electron matter (npe matter) [15j is shown. 

At low temperatures T < O.lkeV the shear viscosity of 
quark matter starts to be significant. Above that tem- 
perature the bulk viscosity is responsible for the stability 
of the star. At temperatures above 1 — lOMeV, where 
the bulk viscosity starts dropping, the region is unstable 
towards GR emission from r-modes. 

The observed distribution of Low Mass X-ray Binaries 
(LMXB's), presumed to be old pulsars spun up by ac- 
cretion of matter from the binary companion to become 
rapid millisecond pulsars, is indicated in the plot by a 
shaded box. It can easily be seen that ordinary neutron 
star models place these pulsars within an unstable region, 
leading us to exclude the classical model of neutron stars. 
On the other hand, the quark matter model places them 
within a stable region. One cannot, however, exclude the 
possibility of npe matter in these stars if one considers 
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T [MeV] 



FIG. 2: (color online) The r-mode instability window for nor- 
mal quark stars (thick lines) for strange quark masses of 100, 
200, and 300 MeV, and neutron stars (thin dotted line) [l5| . 
The shaded box represents the region in which most LMXB's 
are observed. 



a more sophisticated neutron star model with a liquid 
core surrounded by a nearly static crust. This configura- 
tion leads to damping due to viscosity in the boundary 
layer between the oscillating fluid and the crust, which is 
about 10 5 times bigger than the damping from the shear 
viscosity in the interior. This would place the LMXB's 
in a stable region, see, for instance, Ref. [171 ] . 




B. R-mode Instability Window for the CFL Phase 

Color superconductivity changes the situation from 
that of unpaired quark matter. The severity of this 
change depends on the type of pairing and the value of 
the gap. The CFL phase gaps all quark modes. Here we 
consider the contribution of quarks on the bulk viscosity 
and r-mode instabilities of the CFL phase. 

The critical temperature might be as large as T c ~ 
50 MeV when the quark chemical potential /i s ~ 400 
MeV |5Q1 leading to A% FL ~ 89 MeV. The energy 
gap has a dependence on the temperature as Aqfl = 



acfl 



1 



(*)' 



We neglect the effect of color superconductivity on the 
thermodynamic coefficients A i: B i: and Cj. Therefore, 
the effects of the CFL phase are simply to reduce the 
Xi's as in Eq. (fT8|) . Non-leptonic weak interactions have 
two quarks on both the forward and backward channels, 
each of which will be suppressed exponentially, leading 

to H x {A C fl) = H^Acfl) * e- 2 -^. Thus, 



A 1 (Acir i )/Ai(0) = e- 



2& CFL 



(32) 



The two other neutrino-emitting interactions both have 
one quark flavor in both channels. Hence, the reduction 



FIG. 3: (color online) The quark contribution to the bulk 
viscosity as a function of temperature at constant frequency 
"1000 Hz" and r-mode instabilities of the CFL phase for dif- 
ferent values of the strange quark mass. The thin lines show 
the results for unpaired quark matter. 



to the interaction rates will take the form if 2 , 3(^cfl) = 
H2, 3{Acfl) — e C <F L - , Then the reduction of A2, 3 is 



A 2 ,3(Acfl)/A 2 ,3(0) = 



(33) 



We follow the same procedure as for normal quark mat- 
ter to calculate the bulk viscosity and r-mode instability 
window. Since the shear viscosity only participates in 
the r-mode instability at low temperature and all quark 
flavors are strongly suppressed, we can safely neglect the 
quark shear viscosity. The results for a critical tempera- 
ture T c = 50 MeV are shown in Fig. [3] 

Quarks allow for a narrow stability window at high 
temperatures T > 5 MeV. This is due to the fact that 
color superconductivity may actually increase the bulk 
viscosity as can be seen in Fig. [3l 

It is clear from Fig.[3]that considering only the quarks' 
contribution to the stability window places all LMXB in 



an unstable region. However, one cannot yet draw any 
significant astrophysical conclusions; there are other dis- 
sipative phenomena besides just the bulk and shear vis- 
cosities due to quarks in the CFL phase, such as the bulk 
and shear viscosities due to superfluid phonons [Hll |52| , 
and the bulk viscosity due to kaons [53(. Without a 
complete analysis of all these phenomena, no significant 
conclusions can be drawn. However, we do notice that 
at large temperatures there is a narrow stability region. 
Incidentally, this is the region where young pulsars are 
born, which means that young quark stars in the CFL 
phase are produced in a stable region and then have to 
cool down before entering the instability region and slow- 
ing down. 

The result for the r-mode instability window differs sig- 
nificantly from that in Ref . [TtJ , which arises from the fact 
that the bulk viscosity may be enhanced rather than sup- 
pressed due to the suppression of the interaction rates. 
This enhancement was not considered in Ref. 11711. 



C. R-mode Instability Window for the 2SC Phase 

The strange quark mass m s ~ 100 MeV is larger than 
the up and down quark masses m Ut d ^ 10 MeV. There- 
fore, there will be a Fermi surface mismatch between the 
strange quark and the up and down quarks. In this case 
the up and down quarks can pair and form a condensate, 
the so-called 2-flavor color-superconducting (2SC) phase. 

In the 2SC phase, the red and green quarks acquire a 
gap A 2 sc, whereas there is no gap for the blue quark (the 
choice of color is arbitrary) . The critical temperature is 
T c ~ 30 MeV, leading to Ajf c ~ 52.6 MeV [El]. 

Since no global symmetries are broken, the 2SC phase 
is not a superfluid. And since there are ungapped quarks, 
they will dominate all interaction rates and the physical 
properties of the system. They will only be different from 
normal quark matter by a factor as will be shown below, 
see also Refs. (HIH. 

The effect of gap the 2SC phase on the bulk viscosity 
is less severe than in the CFL phase. The reason is that 
there is an ungapped color for the up and down quarks 
and that the strange quarks remain unpaired. 

Since weak interactions do not change the quark colors 
at the vertices, we can decompose the non-leptonic weak 
interaction (the first interaction in Fig[T]) into four differ- 
ent interactions according to the quark colors conserving 
the color at each vertex. This translates to a decomposi- 
tion according to the gapped and ungapped modes of the 
quarks playing a role in the interaction as in Ref. (2o| , see 
Fig. HI 

This leads to a reduction of the interaction rates T\ 
and fi, 

rI ,. . 4 2A 2SC 4 A 2SC 1 

#i(A 2S c) = ge ^ + ge ~ + g> ( 34 ) 

/ a \ 6 A 2SC 3 r~-\ 

#i(A 2SC ) = g e ~ T + 9' ( 35 ) 
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FIG. 4: Contributions to the process u + s-t-^u + din the 
2SC phase. A gapped fermion is marked with the gap A2SC 
at the respective line [25j . 



which gives 



. /a s /> 2 2A 2SC 5 A 2SC 2 ,„„, 

Ai(A 2SC )/A 1 (0) = - e T^ + - e -^ + -. (36) 

This gives a different fraction for the ungapped mode 
compared to that in Ref. [25j. This is due to their as- 
sumption that Sfii > which translated into saying that 
X 1 (A 2 SC) = X 1 (0)H 1 (A 2 SC). Assuming that 5fx t is pos- 
itive neglects the difference between the suppressions of 
both directions of the interaction. 

The same thing can be said about the other two weak 
interactions. However, since there is only one quark par- 
ticipating in the interaction in each channel, we decom- 
pose each interaction into two. the terms A 2j 3 then take 
the form 

A 2 (A 25C )/A 2 (0) = ^ + i e -^, (37) 
A 3 (A 25 c)/A 3 (0) - i + He-T. (38) 

The gapped quarks' contribution to the shear viscosity 
can be neglected, then the reduction to the shear viscos- 
ity is only (5/9) 1 / 3 as in Ref [13 • 

Choosing the critical temperature to be T c = 30 MeV, 
the bulk viscosity and r-mode instability window can 
both be calculated. The results are shown in Fig. 

The results for the r-modes are also different from those 
in Ref. [171 ]. This is due to the assumption that the un- 
gapped mode leads to a 1/9 reduction of the interaction 
rate whereas we have shown that it should be 2/9, and 
due to the fact that the reduction of the interaction rate 
has an effect on the bulk viscosity that is not as trivial 
as only suppressing it. 

It can be seen from Fig. that the effect of the 2SC 
phase is mainly to shift the bulk viscosity and r-mode 
instability window sideways by a multiplicative factor. 
LMXB data do not exclude the possibility of having a 
2SC phase in pulsars. 
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FIG. 5: (Color online) Same as Fig. [3j but for the 2SC phase. 



D. R-mode Instability Window for the CSL Phase 



If the strange quark mass is large enough, the condition 
of charge neutrality requires that the electron chemical 
potential to be large. The strange quark mass and the 
large electron chemical potential will cause a substantial 
difference in particle species' Fermi momenta and pairing 
between different quark flavors may not be possible. The 
remaining option, then, is Cooper pairing of each flavor 
with itself. 

To maintain the fermionic antisymmetry of the Cooper 
pair wave function, single-flavor pairing phases have to be 
either symmetric in color, which greatly weakens or even 
eliminates attractive interactions, or symmetric in spin, 
which causes particles of the same spin to pair together 
and, therefore, the Cooper pairs will have a total spin of 
1. 

For spin-1 color-superconducting phases which are dis- 
cussed in Refs. [E^, [HH, |57[, the situation is very similar 
to the 2SC phase since there are ungapped modes. How- 
ever, it has been suggested [58| that for a particular spin- 
1 phase, the CSL phase, all modes are gapped. There 
are two modes gapped with an energy gap Acsl and 
the "ungapped" mode is now gapped with XiAcsL = 
Jg* Acsl where rrii and fa are the mass, and chemical 
potential of the quark flavor. 

Following the same procedure as for the 2SC phase, 
the weak interaction rates Xi take the form 



Ai(A CS l)/Ai(0) = -e ~ + g e ~ T +g( e T + e T j 

+ — [e t +e t j, (39) 

A 2 (A C5i )/A 2 (0) = i (e-*^+e-*^) + \e-^, (40) 

. /A Wi / r\\ 1 / X d & CSL X » A CSL \ 2 A CSL 

A 3 (A CSL )/A 3 (0) = -(e •? +e r J + - e t . (41) 



Choosing a critical temperature of 5 MeV and ^ = 0.8 
[11], we calculate the bulk viscosity and the critical fre- 
quencies of the CSL phases. The results are shown in 
Fig. [6l The effect of the quark shear viscosity has again 
been neglected since it only acts at low temperature and 
is then strongly suppressed. 

In Fig. [H] we see that the quark contribution places 
LMXB's partially within the window of stability. The 
situation will change, however, if we consider a smaller 
critical temperature than 5MeV. At a small enough T c 
LMXB's would be completely placed in a stable region. 



I 

It must be noted that the CSL phase is a superfluid, 
and so, much like in the CFL phase, there are many 
other dissipative phenomena that need to be considered, 
such as the bulk and shear viscosities due to superfluid 
phonons, bulk viscosity due to Goldstone bosons, and 
friction between both superfluid and normal fluid com- 
ponents (59| . 
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sidered. 

The r-mode instability window for the CSL phase pro- 
vides a wider stability region which ends at a tempera- 
ture T ~ O.OlMeV, see Fig ©. The situation is not fully 
resolved in the case of the CSL phase due to the uncer- 
tainties in the values of the energy gap and the critical 
temperature, and the fact that there are other dissipative 
phenomena to be considered. 

The 2SC phase has a structure very similar to un- 
paired quark matter, see Fig. (O, and the conclusions 
for the unpaired phase can be similarly drawn for the 
2SC phase. This is because in the 2SC phase there is 
always an ungapped quark color. This means that h of 
the quarks behave exactly as the unpaired quark mat- 
ter phase. Moreover, the strange quark always remains 
gapless. This leads to a suppression for the non-leptonic 
weak interaction close to |, for the strange quark Urea 
processes it is | , and for the down quark Urea processes it 
is |. Ref. (l7j suggests that the 2SC phase is marginally 
inconsistent with LMXB's data, we have shown that the 
2SC is completely consistent with the LMXB's observa- 
tions, see Fig. 

To better understand the reasone for the enhancement 
of the bulk viscosity (and of the r-mode stability window) 
due to color superconductivity, we need to look deeper 
into the behavior of the bulk viscosity with respect to the 
oscillation frequency w and the difference in interaction 
rates \. It is hard to see that from Eq. (fT2"|) . Let us 
consider the case of a single-interaction. If only ruled 
by the non-leptonic weak interaction, the bulk viscosity 
takes the form 1411 



FIG. 6: (Color online) Same as Fig. [3] but for the CSL phase. 



VI. DISCUSSION 

In this article, the contribution of quarks to the r-mode 
instabilities was studied. The main conclusion is that one 
cannot neglect the quark contribution to r-mode insta- 
bilities even if all quarks were gapped. 

It is shown in Fig. ([2]) that the r-mode instability win- 
dow is very different for quark matter than for neutron 
star matter. This might provide us with an opportu- 
nity to distinguish between quark stars and neutron stars. 
Further investigation is clearly needed. 

Each different color-superconducting phase has a dif- 
ferent structure for the r-mode instability window. 
Quarks in the CFL phase generate a narrow window of 
stability at temperatures which are, incidentally, the tem- 
peratures at which a young pulsar is formed, see Fig. (|3j). 
This implies that a newly born quark star in the CFL 
phase is generated in a stable region, and has to cool 
down first to reach the instability region which starts at 
T ~ lOMeV. One cannot draw more conclusions because 
the CFL phase has several other contributions to be con- 



Cnon — 2 ' (42) 

•J 2 + (XxAx/n) 

As can be seen from this formula, there are two re- 
gions with two different dependences on the interaction 
rate: 1) The region where uj 2 (\\Ai/n) 2 , in this region 
the bulk viscosity can be approximated by Ci — T^^i 
and any reduction to Ai leads to a similar reduction of 

the bulk viscosity (e f in the CFL phase). 2) The 
region where uj 2 <C (XiAi/n) 2 , in this region the bulk 

nC 2 

viscosity is approximated by £i ~ ^a^ 1 and any reduc- 

tion to the interaction rate (e r in the CFL phase) 
leads to an enhancement of the bulk viscosity (e t ). 

Another way of describing the argument above is to 
say the bulk viscosity is a resonance effect between the 
density oscillation on one side, and the interaction rates 
that try to restore /3-equilibrium during these oscillations 
on the other |25j |. 

The thorough reader might notice that our results for 
the normal quark matter are a bit different -albeit close- 
to those presented in Ref. [ItJ , there are several reasons 
for this discrepancy, unlike Ref. [lT} , which considers the 
bulk viscosity only due to the non-leptonic weak interac- 
tion, the bulk viscosity coefficient considered here takes 
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into account the contribution of both the non-leptonic 
and the Urea processes. Another source of discrepancy 
was due to the assumptions used to calculate the chemical 
potentials of quarks, in here we determine the chemical 
potential of each quark species from the charge neutral- 
ity and /3-equilibrium conditions in a fashion similar to 
Ref. gj. 

This manuscript represents a first step towards a more 
complete treatment of the r-mode instabilities of color- 
superconducting quark matter. As a follow up to this 
work, one could attempt to calculate other dissipative 
phenomena in deconfined quark matter, such as the shear 
viscosity in superfluid phases like the CFL [54] and the 
CSL phase, the bulk viscosity due to superfluidity in the 
CFL [5l( and the CSL phase, see also Ref. [fj(J, or dissi- 
pation due to electron-electron scattering in the system, 
or due to surface friction between a rotating core and a 
nearly static crust in the case of hybrid stars. Supple- 
menting these results with the ones in this article should 
provide us with more accurate figures for the r-mode in- 
stability window in the superfluid phases, which would 
then allow us to draw more firm astrophysical conclu- 
sions. 

In this work a constant density approximation of the 
star was assumed, however, one can, in principle, general- 
ize the calculations to include different models that might 
result in a density profile of the star, which will lead to 
a radis-dependent viscosities and dissipative timescales. 

While preparing this manuscript, an article which ad- 



dresses the dissipation of r-mode instabilities in the CFL 
phase due to the viscosities calculated in Refs. [H2, [53| 
was published, Ref. [6l|. Some discrepancies between 
their treatment and the one presented here need to be 
pointed out, first, in order to calculate the dissipative 
timescale associated with bulk viscosity they used the 
approximation that was used in Ref. [19 . |2l| in which 
the Lagrangian perturbation of the fluid is approximated 
by an Eulerian perturbation. This approximation greatly 
overestimates the dissipative timescale of the bulk viscos- 
ity, as was shown in Refs. Another point worth 
mentioning was that the shear viscosity formula that was 
used is different, it has a temperature dependence of T~ 2 
which was derived in [62j, however, it was shown in [47j 
that a more proper treatment leads to the formula we 
have used in here, i.e. Eq. (fT!))) . 
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